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1. Introduction 

In the standard lattice QCD simulations, large volume, chiral, and continuum limits have to 
be approached in order to obtain reliable physical quantities. That is very challenging even with 
the fastest supercomputers currently available. In most of the previous works, the chiral limit was 
sacrificed and the simulations have been done with considerably heavier up and down quarks, with 
non-chirally symmetric Dirac operators, such as the Wilson or KS-type fermions. 

Recently, an alternative approach, which gives a priority to the chiral limit rather than the 
infinite volume, is pursued based on the understanding of the finite volume effect using chiral 
effective theory. That is the lattice simulation in the e-regime. For such a simulation, the exact 
chiral symmetry for the lattice fermion is essential as one treats extremely small quark masses of 
0(1 MeV). We use the overlap Dirac operator ^ 

D = mo(l + 75Sgn//iy), (1.1) 

which has the exact chiral symmetry through the Ginsparg- Wilson relation Hw = 75Dvi/(— mo) 
is the standard hermitian Wilson-Dirac operator with large negative mass —mo. 
If the size of the box, L, is small and satisfies 

l/AQCD^L<^l/m„, (1.2) 

where Aqcd the QCD scale and niji the pion mass, an expansion in terms of £^ = m^j/Aj/y (Ayy is 
an UV cutoff, e.g. AkFji) is valid after a careful treatment of the zero-momentum modes 0. In 
such a small physical volume, the e-regime, the low-energy constants, such as the chiral condensate 
Z and the pion decay constant Fji, can be extracted through the current correlators. Since the quark 
mass is already very small, the chiral extrapolation is not necessary. 

Most previous lattice simulations in the e-regime were hmited to the quenched approximation 
[§] except for a few pioneering works with rather heavier quarks The new project by the 
JLQCD collaboration started in March 2006 is aiming at performing large scale simulations of 
QCD with the overlap fermion formulations. We also explore the e-regime by pushing the quark 
mass down to a few MeV. We employ topology conserving actions to assure the overlap fermion 
determinant to be smooth [^]. 

For the gauge action we use the Iwasaki action with jS = 2.30 — 2.35. With many additional 
algorithmic efforts we have generated thousands of configurations of two-flavor QCD on a 
16^ X 32 lattice. The lattice spacing is a ~ 0.11 — 0.125 fm. As an exploratory run we attempted 
to reduce the quark mass down to 2 MeV at jS = 2.35, at which L ~ 1.8fm. With this small quark 
mass the condition ( ll.2| ) is safely satisfied. 

In this report we concentrate on the study with the smallest sea quark mass ma = 0.002 (m ~ 
2 MeV) at 2 = 0. We first report on the numerical cost for this simulation in Sec. ^, and describe 
the determination of the lattice spacing in Sec. |3[ Preliminary results for the low-energy constants, 
£ and Fji, are obtained through the low-lying eigenvalues (Sec. Q) and through the pion correlators 
(Sec. The summary and discussion are given in Sec. ^. 
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Figure 1: The history of CG iteration along the HMC steps. The sea quark mass dependence seems weak 
when ma < 0.03, and the numerical cost with ma — 0.002 (red) is almost the same as the case with ma — 0.02 
(blue). 



2. Numerical cost 

Naively, the simulation cost increases for small quark masses as an inverse power of the quark 
mass m, and tends to diverge near the massless limit. It is, however, not the case in the finite 
volume. The numerical cost or the iteration count for the solver, say the conjugate gradient (CG) 
solver, is roughly proportional to the condition number 

\Kax + m\^/\^in+m\^, (2.1) 

where Xmm,max denotes the lowest (highest) eigenvalue of the overlap Dirac-operator. In the large- 
volume limit, the numerical cost is primarily determined by the quark mass since A„„>, ~ 1 I'D/ ~ 
and Xmax is insensitive to the gauge configuration. In the small quark mass limit at a fixed V , on 
the other hand, m is made smaller than 1 /ZV and {Xmin) — z/^V determines the condition number 
(2.1) ^ Here, the numerical factor z can be estimated using the chiral random matrix theory. For 



Nf = 2 and 2 = 0, it is z ~ 4.34. 

In fact, as Fig. [l] shows, the quark mass dependence of the iteration count is weak for am < 
0.03, and the simulation cost for am = 0.002, which corresponds to m ~ 2 MeV (with an assumption 
for the renormalization factor = 1.8 as in the quenched theory) is almost the same as that for 
ten times heavier quark mass am = 0.020. However, we should note that the auto-correlation time 
seems longer for smaller quark mass. We do not have enough statistics to precisely calculate, but 
we expect it is O(IOO) trajectories. 

On a half-rack (512 nodes) of the IBM BlueGene/L (2.8 Tflops peak performance), we need 
roughly one hour per trajectory. In the following analysis, we use 100 configurations in 2 = 
topological sector sampled from 1400 trajectories (first 400 are discarded for the thermalization). 
For the details of our numerical simulations, we refer to the contribution by Matsufuru [P]. 



When there are zero modes, A„„„ = 0, one must expUcitly subtract them from the determinant. 
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Figure 2: The static quark potential (left) and the lattice spacing with different quark mass (right). The 
quark potential has no indication of string breaking and a ^ 0.11 can be obtained, which is also consistent 
with the chiral extrapolation from the heavier quark mass points, as the dotted line in the right panel shows. 



3. Static quark potential with almost massless sea quarks 

To determine the lattice spacing a we measure the static quark potential. In a very small quark 
mass regime one might expect the string breaking, but our results in Fig. ^ (left panel) show no 
indication. This is probably because the volume is too small to contain two static-light mesons. 
The overlap of the Wilson loop with the two static-light states could also be a problem. 

We calculate the Sommer scale ro as usual and obtain a ~ 0.111(2) fm (with an input ro = 
0.49 fm), which is consistent with the chiral extrapolation from the heavier quark mass points as 
shown in the right panel. 



4. Low eigenmodes 

We expect that the eigenvalue distribution of the Dirac operator is well described by the chiral 
random matrix theory (ChRMT) in the £-regime. We compare them with our lattice data. 

The eigenvalues of the overlap-Dirac operator lie on a circle in the complex plane, and we 
project them as 

where X"^' is the i-th complex eigenvalues of the overlap-Dirac operator D. Note that A; is very 
close to ImAf^ for the low-lying modes. 

ChRMT predicts that the lowest mode with an eigenvalue Ai in Nf = 2 theory feels strong 
repulsive force from zero. That is more prominent than the quenched case, and numerically one 



obtains [10] 



(Ai)ry = 1.77 (A^/ = 0), (Ai)ry =4.34 (A^/ = 2), (4.2) 

where the chiral condensate Z could depend on the number of flavors. The HMC history of Ai 
shown in Fig. ^ is consistent with this theoretical expectation. The lowest eigenvalue is pushed up 
more for lighter quarks than for heavier quarks. 
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Figure 4: Ratio of the low-lying eigenvalues. "///' denotes Xi/Xj. Left: we plot the data with {Nf,Q) ~ 
(0,0), (2,0), (0,2) (red crosses). The ChRMT predictions are shown by blue stars. Right: quark mass 
dependence of the ratios "2/1", "3/1" and "4/1". Crosses are ChRMT results with {Nf,Q) = (2,0) = (0,2) on 
the left, and (A^^, Q) — (0,0) on the right, where we also plot the corresponding quenched results (squares). 
Circles are lattice data in two-flavor QCD. 



Since the value of chiral condensate cannot be directly compared for different sea quark 
masses, it is convenient to look at the dimensionless ratio of the eigenvalues Xj/Xj. We present 
some of such ratios of low-lying eigenvalues on the left panel of Fig. In the plot, the results are 
shown for both quenched and unquenched (Nf = 2) cases at 2 = 0, as well as for the quenched data 
at 2 = 2. All the data show a good agreement with ChRMT (blue star symbols). ChRMT predicts 
that the {Nf, Q) = (2,0) theory has the same low-lying eigenvalue spectrum as the {Nf, Q) = (0,2) 
theory, which is well reproduced by the lattice data. 

It is also interesting to see the quark mass dependence of the ratios (the right panel of Fig. 
The data with am > 0.03 are close to the quenched ChRMT, which is shown by the crosses on the 
right end, and below am = 0.03 they suddenly drop to the Nf = 2 results (crosses on the left end). 

We obtain the lattice bare value of (E^/=2) V3 - 228.9(3.6) MeV from Ai measured in am = 
0.002 configurations. Here, the renormalization of £ is not taken into account. 
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Figure 5: The left panel shows the effect of low-mode averaging (LMA) on the pion propagator(for one 
configuration): data with (blue stars) and without (red crosses) LMA are shown for comparison. A quadratic 
fit is attempted for the LMA pion propagator on the right panel. 



5. Pion correlators in the £-regime 



The pion correlators in the £-regime are largely deformed by the finite volume effects, and no 
longer an exponential but a quadratic function of time t, 

J d\{P"{x,t)P\0,0)f = 5"'''[A{t/T -\/2f+B], 



2L 



Ml 



,B 



mT 



(5.1) 
(5.2) 



where pL = niLV and ji' = /i(l + 3j3i/2F^V'/^). The effective chiral condensate 'LQ{ii) has strong 
dependences on Q and m, which are expressed by the modified Bessel functions (see [Q] for the 
details). Here, (0)^ denotes the expectation value in the Q topological sector. 

We calculate the pion correlators on our lattice using the low-mode averaging (LMA) tech- 



nique [11 1. We find that ~90% contribution comes from the lowest 100 modes of the Dirac opera- 
tor, and the LMA is very effective, as the left panel of Fig. ^ shows. 

With an input 1)1^ = 229(4) MeV from Ai, the quadratic fit to the correlator works well as 
shown on the right panel. We obtain a preliminary result F,^ = 86(7) MeV with a fitting range 
t = [10, 22], for which ;t;7dof ~ 0.25. The resuh is very preliminary, not only because the statistics 
is hmited so far but detailed analysis of the systematic errors is yet to be performed. 



6. Summary 

We pushed the dynamical overlap fermion simulation towards the chiral limit and in fact suc- 
ceeded a simulation at m ~ 2 MeV on the lattice with L ~ 1.8 fm and a ~ 0.11 fm. On the finite 
volume lattice, the CG count does not diverge in the chiral limit as expected from the chiral effective 
theory. The distribution of low-lying eigenvalues of D is consistent with the ChRMT expectation, 
from which we may extract the value of the chiral condensate. 
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The measurement can be extended to the hadron correlators. As a first test, we calculate the 
pion correlator and find that the LMA works extremely well in the e-regime. Precise determination 
of Fj[ is feasible with such calculations. 

So far, our simulation has been done at a single parameter set, Q = and am = 0.002. To study 
other topological sectors and different quark masses is very important, because in the £-regime Q 
and m dependences are prominent and help to improve the sensitivity to the low-energy constants. 

Numerical studies are performed on Hitachi SRI 1000 and IBM Blue Gene at High Energy Ac- 
celerator Research Organization (KEK) under a support of its Large Scale Simulation Program (No. 
06-13). Some parts of numerical analysis were carried out on SX8 at YITP in Kyoto University. 
This work is supported in part by the Grant-in-Aid of the Ministry of Education (No. 13135213, 
16740156, 17340066, 17740171, 18034011, 18340075, 18740167). 
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